Heavy Quark Potentials in Some Renormalization Group Revised AdS/QCD Models 
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We construct some AdS/QCD models by the systematic procedure of GKN. These models reflect 
three rather different asymptotics the gauge theory beta functions approach at the infrared region, 
/3 oc — A^, — and P oc —A, where A is the 't Hooft coupling constant. We then calculate the heavy 
quark potentials in these models by holographic methods and find that they can more consistently 
fit the lattice data relative to the usual models which do not include the renormalization group 
improving effects. But only use the lattice QCD heavy quark potentials as constrains, we cannot 
distinguish which kind of infrared asymptotics is the better one. 

PACS numbers: ll.lO.Hi, 11.25.Tq, 12.38.Aw, 12.38.Lg 



I. INTRODUCTION 

The Duality between string theories in Anti-deSitter 
space and conformal field theories on its boundary, i.e., 
AdS/CFT, are powerful tools for understanding the 
strong coupling gauge theory phenomena. In its most 
well known formulations, the duality is between the 4- 
dimensional A/" = 4 super Yang-Mills field theory and 
type IIB super string theory in the AdS5 x ^5 background 
[H, In practice, our greatest interest is not A/" = 4 
SYM, which is conformal and non-confining, instead it is 
the Nc = 3, Nf = 4 or Nf = 6 Quantum Chromodynam- 
ics, i.e. QCD, which has running coupling and confining 
properties. 

To obtain predictions about the practical QCD the- 
ories through AdS/CFT, or more generally, the holo- 
graphic principle, people develop two ways, (i) the top- 
down method and (ii) the bottom-up approach. In the 
first method, people either construct dual string/gravity 
descriptions for gauge theories with running couplings, 
e.g. cascading gauge theories^, or build models for 
gauge theories with fewer super-symmetries relative to 
the A/ = 4 SYM[3|0. For examples, some models im- 
plement fundamental flavors by adding probe branes in 
various exact or asymptotically AdS^ x background 
QQ, see also Q for reviews. The top-down method 
preserves fundamental structures of string theory and its 
dual gravity description is 10 dimensional. 

In the second way, now known as AdS/QCD, people 
start from some 5 dimensional effective gravity theory 
and calculate the relevant 4D gauge theory quantities 
through holographic method. By requiring the results 
coincide with those from QCD phenomenologies, people 
obtain the general properties of the dual gravity back- 
ground, which can be refined for further phenomenolog- 
ical goals. In the earliest iniplementation of this ideal 
by Polchinsky and Strassler [9|, an infra-red cut off on 
the dual AdS^, background is introduced to implement 
confinement, hence the name hard wall model. Since 



the space-time of gravity descriptions ends below the 
cut off, appropriate boundary conditions must be im- 
posed on this point by hand, see e.g., reference [l3] and 
[ll| . In reference [T^l, it is proposed that, by adding 
quark masses/condensates, this model can become a re- 
markably good model of chiral symmetry broken dynam- 
ics. Reference [l^ considered a natural extension of this 
model by replacing the pure AdS background with a 
dilaton-deformed one and get very similar results. To 
improve the Regge behavior of highly excited rho mesons 
and higher spin mesons in the hard wall model, reference 
proposed the soft wall model, where the infra-red 
hard wall is replaced by a dilaton quadratically depend- 
ing on the holographic radial coordinate. While refer- 
ence [3 and [Tg proposed to include the same effects 
by a simple factor e'^^ in front of the dual geometry's 
metric. Different from this guess and trial approach, in 
references [l3|, [Hi and [11] GKN(Giirsoy, Kiritsis, and 
Nitti) propose a systematic procedure for dual geometry's 
construction and provide several models as illustrations. 

We will introduce GKN's procedure and construct our 
own models in the next section. To make our dual geome- 
tries as simple as possible, so that more QCD quantities 
can be calculated conveniently and analytically, we do 
not require our models' /3 function to exactly coincide 
with the perturbative QCD theory to higher loops, like 
those provided in GKN's original work. The next next 
section discusses the calculation of heavy quark poten- 
tials in the resulting models. Numerical fittings as well 
as discussions about the quantum corrections are also 
provided in that section. The last section contains our 
main conclusions and prospects for future studies. 



II. THE MODELS 

By GKN's arguments, the non-critical string back- 
ground dual to QCD like gauge theories can be described 
by the following Einstein frame action 
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where we have neglected the effects of an axion field, 
which is of 0{N~^) relative to the terms written out 
explicitly; the potential V{^) encodes various contribu- 
tions such as the higher a' corrections and the integra- 
tion of RR 4-form field whose flux seeds the D3 branes 
and U (Nc) gauge group. Its functional form cannot be 
derived out explicitly from the first principle, but can 
be determined through inputs of QCD phenomenologics, 
e.g., the P function describing the running of coupling 
constant. 

Setting the dual gravity geometry as 



(i){u) (2a) 
e2'4(-dt2 + dx ■ dx) + du^, (2b) 



and introducing phase space variable X = , GKN find 
that the dilaton field's equation of motion and Einstein 
equations following from the action ([1]) can be written as 
first order forms. 




!/($) = Vo{l-X-}e 
^ and A' = ^. 

du du 



Xd<S> 



(3a) 

(3b) 
(3c) 



where = ^ and A' = After identifying the 

exponentiated dilaton field with the 't Hooft coupling 
constant A = NcQym and the metric function e"**-"^ with 
the QCD energy. 



A 



cx E 



(4) 



GKN set up an explicit relation between the gauge theory 
(3 function and the gravity theory quantities $ and A, 



(5) 



To implement asymptotical freedom in the dual gauge 
theory, it is required that A{u) — > it, <I> — > — oo as 
M — > oo, i.e., the dual space-time is asymptotically anti- 
deSitter. Obviously, in this construction of AdS/QCD 
models, the non-perturbative (3 function determines the 
whole structure of the dual gravity theories. In turn, 
as long as we know the latter, we can use it to calcu- 
late more QCD phenomenological quantities through the 
holographic method. 

In this paper, we will consider 3 models, each with 
exact (non-pcrturbativc effects included) (3 function 



Model 1 
Model 2 

Model 3 
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-bX^ (6) 
-boX^ - ^^lA^ , 6o , 6i > (7) 
boX' 



1 + biX 



bo , bi>0 



(8) 



respectively. At the perturbative limit, model 1 repro- 
duces QCD's beta function to 1 loop order; model 2 re- 
produces the desired (3 function to 2 loop order; model 3 



only reproduces the desired(perturbative) {3 funtion to 1 
loop order. Nevertheless, there are evidences that non- 
perturbative effects indeed produce f3 functions of the 
model 3 type^see e.g. ^ ,[2l| and The author 

of reference [23 tells us that, this kind of beta function 
reproduces exactly all of the known results and also the 
two loops. At the strong coupling limit, model I's /3 func- 
tion approaches infinite as — A^, model 2 as — A"^, while 
model 3, —A. By GKN's general confinement criteria, 
model 3 is critical confining, while model 1 and model 2 
are super-confining. 

For model 1, according to GKN's procedure. 
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6A" = — ^ A= — 



dA 



bX 
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-r^bX^e^''^ (10) 



(11) 



so the dual gravity background can be written out ex- 
plicitly 



ds^ = e^'^i-dr + dx ■ dx) + te^'^dX^ (12) 
$(A) = logA, e^ = ei^ , e^ = -^ (13) 

where we use A as the radial holographic coordinate. Us- 
ing eq (|3c)) . we know that the potential of this model's 
dilaton field is, 



V{X) = l/o(l 



6^ 
9 



)e 



(14) 



From this potential, we easily see that the dual gravity 
background is asymptotically AdS in the ultraviolet limit 
but not so in the infrared region. 
From equation (fTTj) . we know that 



du 



6A2 



-dA 



(15) 



By definition and iteration formulas of the incomplete 
Gamma function, this leads to 



oo g-|6A 



ilbxy 



dihx)^^e-T[-iAbx\ 



'l-(l6A)-Vr[0,l5A] 



(16) 



Note that xr[0, x].j;^o 0, relative to the first term, the 
second term on the r.h.s of (|16p can be neglected. So in 
the ultraviolet region. A— s-0,u— s-^,A=jj^j. Us- 
ing domain wall coordinates (|2bp , the gravity background 
can be easily looked out being asymptotically AdS. While 
in the infra-red region A — > cx), the scale factor of the do- 
main wall coordinate = — > 1 this forms a natural 
softwall for the dual gravity background. 
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For model 2, by the same procedure as previous we 
know the dual gravity background and the corresponding 
dilaton potential are 



dl 



e^''{~dt^ + dx-dx)+te 
bo 



2„2D 



log A 



6oA2 + 5iA3 



dX^ (17) 



(18) 



and 

V{X) = K)[l-^(^oA + 6iA2)2]ef(^o^+l''i^') (19) 

respectively. Obviously, in the ultraviolet region A ^ 0, 
this model has the same asymptotical background — AdS 
space-time — as model 1. While in the infra-red region, 



du 
dX 



6iA3 



2 n 



e-i''^^^(^6iA2)-Vr[0,^6iA2]](20) 



So in this model, the rate of u — > as A — > oo is differ- 
ent from that in model 1 , but the qualitative conclusions 
that, (i) u as A — > GO, (ii) the dual gravity back- 
ground is not asymptotically AdS in the infra-red limit, 
(iii) the dual gravity space-time's scale factor e^_^^ 1 
are common between the two models. 
For model 3, 



1 + biX 



dX 
dA 



^ds^ ^ e^^{-dt^ + dx-dx)+fe^°d}? (21) 

1 6l 

ci>(A) 



log A 



— e''o^ (boX) ''0 
(I + 61A)' W 



V{X) = Vo[l 



blX^ 



9(l + 6iA) 



](l + 6iA)W 



(22) 
(23) 
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In the ultraviolet region, A ^ 0, 

u — > e^. So, like the previous two models, 
this model also has asymptotically AdS geometry. But 
in the infra-red region, this model's dual space-time scale 

factor e (boX) ''o shrinks to zero as A ^ 00, the 
corresponding domain wall coordinate u approaches zero 
as A —> 00 by the following law 



(24) 



In summary, all three models have the same ultraviolet 
asymptotic geometries — AdS space-time plus running 



dilaton field. In the infrared region, the first two model's 
dual space-time scale factor approaches to 1 as A ^ 00, 
but third model's scale factor shrinks to zero size as A 
00. The domain wall coordinate u of model 1 approaches 
zero at the rate u ~ 2e~^^^, that of model 2 at the rate 

u ~ 4e"i''l^^ while that the model 3, u ~ (61 A)"^. 



III. HEAVY QUARK POTENTIALS 

The heavy quark potential is an important quantity 
related with confinement. It is measured with high pre- 
cision in lattice simulations, see reference [2^ and [2^ . 
The results are usually fit into rational polynomials 



E{p) = — + — + ^ 



C 



(25) 



where p is the distance between the quark and anti-quark 
while K, k' and a? are parameters to be determined by 
fittings, C is an irrelevant normalizing constant. In gauge 
theories, this potential is related to the expectation value 
of rectangular Wilson loops of width p and length T 



{W) = eM-E{p)T] 



(26) 



here T can be understood as the time the quark and anti- 
quark are bounded together. By the gauge/string duality 
conjectures, {W) equates to the partition function of an 
open string which is living on the AdS background and 
whose world sheet boundary coincides with the Wilson 
loop 



{W) 



exp[-£;(p)r] ^ Z = j VXf" cxp(-S') 



(27) 



where S is the Nambu-Goto action of the string. Let 
x^ {a, t) denote the classical string profile and ^'^ the 
quantum fluctuation around it. Obviously, the path in- 
tegral will be dominated by the classical string configu- 
ration since it minimizes S, this means 



Z = cxp(-5c[x]) J V£_exp[-S{x 



+ + Sc{x)] 



CXp(-5eN) / V^CXpi^C^OaC] 



(28) 



So in the classical approximation, the heavy quark po- 
tentials read 



(29) 



(30) 



E{p)T^Sc[x] 
While to first order quantum corrections 

E{p)T ^ Sc[x] ^J2^ogdctOa 

In the following, we will calculate the heavy quark poten- 
tials quantitatively in the classical approximations and 
study the quantum corrections qualitatively. We do not 
get quantitative results on the latter. But even in the 
classical approximation level, the models provided in this 
paper predict heavy quark potentials more consistently 
fitting with the lattice QCD results than almost all the 
models considered in 12811. 
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A. Classical Approximations 

In the classical approximations, to have the given Wil- 
son loop as the world sheet boundary, the string must 
have its two ends fixed on the boundary and the center 
part dips into the bulks of the background AdS space- 
time 



, x"^ = o- , A A((t) , x"^'^ fixed 
A(-f) = A(^) = 0, A(0) = Ao 



(31a) 
(31b) 

(31c) 



For this configuration, the classical Nambu-Goto action 
can be written a,a Sc = ■ J dr, with 



1 

na' 



da 







where A^o- denotes 
minimizes Sr we know 



8* 4A 

e 3 e 



g^e2A(g2A+£2g2_D^2^) (32) 

Since the classical configuration 



const. 



e 3 e 



2A(g2A_^^2g2Z?;^2^) 



(33) 



±e 



A-D 



g3<E'g4A 



1 (34) 



in the above derivation, we have used the fact that at the 
center point of the string ct = 0, A^o- = 0. Using eqs ([33l) 
and (p4|) , both the distance p between the two ends of the 
string and its energy Ec can be expressed as functions of 
Ao. 



Er 



2£e- 



dX 



,D-3A^- 



dX 



Vi- 


A-fe-4^ 








- A-ie-4^ 



(35) 
(36) 



A S [0, 1], the parameter c must be less than |. Through 
qualitative analysis of, we know that p is a monotonically 
increasing function of c. 

Under the long separation limit which corresponds to 
c ^ I , the potentials can be calculated easily 



E, 



c,gq,/5— >oo 



2TTa' 



\3p2Ao 



dAA^e^+^ 



eiire^^A 4c. 



27ra' '3' 9 

While at the short separation limit c 0. 



2£e-- (0.596 -I- 0.0471c + 0[c^]) 



Er_ 



(-0.602 0.042c -hO[c2 



(39) 
(40) 



Which means that 
e{c) = 



2f^ t{c%=c(p) 
a'&4/3 p 

ro.359-h 0.003c- 



(41) 



the functional form of c(p) should be obtained by revers- 
ing the relation (|39|) . For example, to first order ap- 
proximation in [log^] , we can use iteration method 
to obtain 



c « ( log [ 



1.19^ 0.0942^,, 1.19£,-1t\-i 
+ (log ) ] (42) 



Since the series expression for p(c) in ea (j39[) is incom- 
plete, its inverse (|42|) is not well-defined as p > 1.19£. 
We assume that c(p)p>o.8/Gcv = c(0.8/Gev). Combining 
eas ([55)) and (|1T|) together, the classically approximated 
heavy quark potential of model 1 can be written as fol- 
lows 



where A = A — Aq while A = Since the string is 
infinitely long, the total energy in ([55]) is divergent. In 
the Gauge theory side, this means that quark-antiquark 
pairs are infinitely heavy. To obtain the finite interaction 
energy, we will subtract from E^. the infinite rest mass 
part, which corresponds to that of two straight strings, 
see reference and [1^ . Taking model 1 as an example, 
this means 



Er 



lX^e^'> 



dX- 



A3 e 



D+A 



1 



Vl-A-fe-4(A-i-i)/c 



dAA^e^+^ 



(37) 



where A = -f-. C:=6An and the 



' still means an ex- 
pression the same as the denominator. To assure the inte- 
grand of the right hand side take real values in the range 



2£2 e(c),=,(p) 



'63 



'&4/3 



27ra'53 



ill^r[i -] 

ira'b's 3 3 



(43) 



Figure [T] compares this result with the direct numerical 
integrations intuitively. From the figure we see that the 
analytical formula captures the main properties of the 
integration remarkably well. For model 2 and model 3, 
the potentials can be calculated similarly, and the results 
are also similar. But the precision of the analytical ap- 
proximation is lower than model 1. For this reason, in 
the following fittings, we use numerics to get the p^Ec.qq 
relations for all the three models directly. 

Ea ipS]) does not have the the standard Cornell po- 
tential form, for its ^ term(or the Luscher term) 
has a p-dependent coefficient, which is proportional to 



5 
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FIG. 1: Heavy quark potentials in model 1, the continuous 
line is the result of analytical approximation (|43|l . while the 
dashed line comes from direct numerical integration of ea (|35p 
and 133. 



[lof 



1.19^1 



This is determined by asymptotical behav- 



ior ^ of the beta function /? cx — and is a common pre- 
diction of all the renormalization group revised models 
constructed by GKN's general procedure [ll], Just 
as we will show in the following, this feature makes these 
models predict p — Eqq relations more consistently fitting 
with the lattice data than those which do not include the 
renormalization effects. 



B. Comparing with lattice results 

Figure [2] compares the the heavy quark potentials fol- 
lowing from our AdS/QCD models and those of lattice 
results. The relevant parameters are listed in table |l] and 
are fixed by fits, which is required to minimize the fol- 
lowing variable 



E 

1=1 



^^AdS/QCD _ ^latt.Y 



(46) 



where the sum is over all the used data points and cr,; is 
the absolute error of the ith point. From the figure we 
see that almost all the three models fit with the lattice 
data equally well. This means that in the AdS/QCD 
framework, only by fitting with lattice data on heavy 
quark potentials, we cannot distinct the infrared behavior 
of the gauge theory (3 function. 



^ If we assume that 

/3 -> ~bXP , p > 1 (44) 
then, by the same calculations as above, we will find that 



[logi^] 3(p-l) 



+ other terms 



(45) 





Pcut 


parameters 


Best fit 


Model 1 


1.0 


4 


4.99 






£ 


4.71 






const 


-0.013 


Model 2 


1.0 


a'bl 


4.38 






e 


4.37 






hi /hi 


0.011 






const 


0.018 


Model 3 


1.0 


a'hl 


4.55 






I 


4.25 






hi /ho 


0.058 






const 


-0.029 



TABLE I: Best fit values for the parameters involved in our 
AdS/QCD models. All dimensional parameters use Gev or 



Gev as units. 
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FIG. 2: Potentials Eqq{p) predicted by our AdS/QCD models. 
The discrete data points comes from the lattice results of 
quenched QCD, reference ^31j]. Since the three models almost 
fit the lattice data equally well, the three lines coincide. 



Taking our FIG. [2] and the Fig.l of reference [111 as a 
comparison, we easily see that our models can be more 
consistently fitted with the lattice data than almost all 
the previous models which do not include the renormal- 
ization group improvements. To make this point more 



clearly, we display the x l^d.o.j 



Pcut 



in figure [3l Our 



Pcut has the same meaning as that of reference [28| , which 
is the length scale below which the lattice data is not used 
to make fits. This is a manufacturing selection. Our mod- 
els can give very low l^d.o.f even when the pcut is set 
as low as lGev~^, while most of the AdS/QCD models 
which do not include the renormalization group effects 
can only do so when pcut is set as 4 — 5Gev~^. In a an- 
other word, to make a good fit, the usual models must 
cut more than one half of the lattice data points while 
our models almost need no cut at all. 

This is a really exciting result and it can be attributed 
to the fact that our models predict heavy quark poten- 
tials which at the short separation limit behaves like 



Ea. 



3 /p. Considering the fact that al- 



6 




Pcut/Ge 



FIG. 3: Variation of the minimal /Nd.o.f with respect to 
pcut- The red, green and blue lines are the results for our 
model 1, model 2 and model 3 respectively. 



most all the AdS/QCD models which do not include the 
renormalization group effects predict Egg p_,Q — > 1/p and 
the quantum corrections do not change it remarkably [26|. 
a natural question arises, will the quantum corrections in 
our models change the short separation hmit of our po- 
tentials remarkably? We change in the following to the 
quantum corrections of the heavy quark potentials. 



C. Quadratic quantum corrections 

Including quantum corrections, to first first order ap- 
proximation, the heavy quark potentials can be read out 
from eqsdUl) and (pO)) . 



(47) 



E{p)T = S,[x] -J2^og {dot Oa) 



To get quantitative results, we need to run 3 steps [26|. 
(i) choose an appropriate gauge and define the fluctu- 
ations properly, (ii) work out the corresponding action 
S[xc + £,] — S[xc] describing the fluctuations and read 
out from it the relevant second order differential oper- 
ators Oa, (iii) calculate eigenvalues of these operators 
and use some regularization procedures to calculate the 
functional determinants of them and finally extract the 
resulting corrections to the potential. 

We will consider the fluctuations in the fixed A — fix it 
on the classical profile — gauge. In this gauge the fluctu- 
ating quantities are ^ and 77, whose boundary properties 
and relations to the complete embedding functions X'^ 
are 

X2^3 = r;2'3(r,a) (48a) 
C(r,a)r=o,T -0, ^{T,a),^o,±p/2 = (48b) 

T]{t, cr)r=0,T = , 7]{t, Cr)<j^o,±p/2 = (48c) 

where X{a) describes the classical string profile which sat- 
isfies cq ((M)) . Substituting these embcddings into eq 



and extract terms quadratic in quantities ^ and 77, we 
have 



1 



(2) 



Gs ( r~is 



/2 ^2 



Gti^xx^'^ + GyyiG^x + Gx\X'^)rf' + Gf^GyyTj'^ (49) 

where the over-dots denote derivatives with respect to 
T and primes denote derivatives with respect to ct. After 
integrating by parts and simplifying by the classical e.o.m 
5(2) can be written in the form ^^O^^+rj^Orjrj, with 



Or, 



4A 



-(A|e^^")[-j9| + (A3e 
-{xle'^")[^dl+~Xh'^d!] 



(50) 
(51) 



where x is the normalized x ~ x ■ p, this means that as 
X varies in the range < x < p, x varies in the range 
Note that 77 is a two component vector 

and similarly 



< 



< 1. 



and 1)^ actually means ?y(2)j^(2) 



77- 



^(2)'^(2)' _^ jj{3)'^{3)'_ 



Reference [26| proved that, for any operator of the form 
0[A, B] ^ B^O^ + A^Ft{v)d^ (52) 



if the parameter A, B encode all the dependences on the 
inter-quark distance p, while v, Ft and are variables 
or operators independent of p, then the functional deter- 
minant of 0[A, B\ must be of the form 

log(det 0[A, B]) = ^ log(dct 1]) + o{T) (53) 

i.e. p dependence of the logdet of such operators is com- 
pletely contained in the coefficient For almost all the 
models considered by [1^, say the typical AdS^/CFT^ 
case, this leads to that the quadratic quantum correction 
to the heavy quark potentials are proportional to p~^. 
However, in the current case, if we write the operators of 
eas ([50|) and ([51]) into the form of eq ([52|) . then the quan- 
tity Ft{v) = A^e^'^ appearing in them will depend on the 
intcrquark distance implicitly, see figure |4] for details. 
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— p = 


O.Oli 


r 




O.W£ 


1 


— p = 


0.21£ 


I'l 






// 


















1 



1 

0.2 0." 



Ft{v) 




0.6 0.1 

V 



1.0 0.2 0.4 0.6 0.8 

V 



1.0 



FIG. 4: Left, the function Ft{v) appearing in the Or, operator 
of eas (|5ip (write it into the form of (|52|l ). the function depends 
on the interquark distance p. Right, Ftiv) for the simple 
AdSzi/CFTi correspondence, the function is independent of 
the interquark distance p 
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A = CXD 



A = CO 



A = n 



www 



A = CO 



(a) 



(b) 



(c) 



FIG. 5: String configurations contributing to the Wilson 
loop correlator. The upper lines denote the boundary of the 
asymptotically AdS space, while the lower dashed lines denote 
its center point. According to reference [23], figure b) is the 
dominant contributor to the Wilson- Polyakov loop correlators 
at the short separation limit; figure c) is the dominator at the 
long separation limit; figure a) displays quantum corrections 
due to the graviton-dilaton exchange process. 



This is both a good and a bad news for us. The 
good news is, it means that we can almost definitely con- 
clude that the quantum corrections to the heavy quark 
potentials is not of the 1/p type. The bad news is, it 
means we cannot derive out explicit functional forms for 
the p-dependence of the quadratic quantum corrections 
through this simple method. To get quantitative results, 
we need to find out the eigenvalue spectrum of and 
On and their functional determinant by numerics. Since 
they cannot be find out analytically, the regularization of 
the resulting functional determinant becomes a big prob- 
lem. At this point, we have no good method to overcome 
this obstacle numerically. 



Reference [27| studied the Wilson-Polyakov loop cor- 
relators in the finite temperature A/" = 4 SYM plasmas, 
from which the relevant heavy quark potentials can be ex- 
tracted conveniently. They pointed out that although at 
the short separation limit, the heavy quark potentials is 
dominated by smooth the string configuration which min- 
imizes the classical Nambu-Goto action, part(c) of FIG|ni 
at the long separation limit, the potentials is dominated 
by string profiles which involves the exchange of the light- 
est graviton-dilaton modes, see the part(c) of FIGO Ac- 
cording to this results, our above considerations of the 
quantum fluctuations (around the classical string profile) 
correction is not the only contribution to the interquark- 
potentials. We should also consider diagrams of the form 
like the part (a) of FIG [5] But this kind of string profile 
contributing to the heavy quark potentials always carries 
a factor of relative to the classical profile does, and to 
make the area of the string world sheet as small as pos- 
sible, the string cannot not fall very deep into the bulk 
of the space-time, i.e. they are mainly around the A — > 
region. So there is no need to worry that this kind of 
correction will change the short separation limit of our 
potentials. 



IV. CONCLUSIONS 



From the previous section, we know that the heavy 
quark potentials following from all our three models 

1.19^4/3. 



asymptote to the same form Eqq oc — l/(p[log■ 



at the short separation limit. This is determined by the 
UV property of the gauge theory beta function /3 cx — A^ 
and is a common feature for all the renornialization group 
improved AdS/QCD models constructed by GKN's pro- 
cedure. Although this is different from the usual Cor- 
nell potential, through numerical fittings and compari- 
son with reference [28|, we know that our models can 
more consistently fit the lattice data, relative to almost 
all the previous AdS/QCD models which do not include 
the renormalization group effects, for example, those of 
[13, [H and 

Numerical results also tells us that, though we can con- 
struct AdS/QCD models with very different infrared (3 
function, for example (3 — A^,— A, — A"^, only by com- 
parison with the lattice QCD heavy quark potentials, 
we cannot distinguish which one is the better. Since at 
long separation limit we neglected the graviton-dilaton 
exchange contributions to the heavy quark potentials, we 
should be more careful on this conclusion. Just as ref- 
erence [igj showed and our numerics indicated, neglect- 
ing the graviton-dilaton exchange contributions does not 
affect the qualitative conclusion that the potentials fol- 
lowing from our AdS/QCD models are linearly confined. 
(Reference [23l pointed out that there is many uncer- 
tainties on the quantitative calculation of the graviton- 
dilaton exchange contributions at the long separation 
limit) Nevertheless, further studies on this point and 
comparison with both the lattice data and the heavy 
quark effective field theories [s^l will be interesting top- 
ics. 

The above content consists the main conclusions of this 
paper. On the model building part, we have no unex- 
pected conclusions. A little summary can be given as 
follows. Three AdS/QCD models are built by the sys- 
tematic procedure of GKN. At the perturbative limit, the 
first model has gauge theory (i function coincide with the 
usual QCD to 1 loop order; model 2's /3 function coin- 
cides with the desiring to 2 loops; while model 3, only 
reproduces the desiring to 1 loop order. Nevertheless, 
due to non-perturbative effects, the (3 function of model 
3 may be the one most closely related to the real QCD 
theory. At the strong coupling limit, model Vs (i func- 
tion approaches infinite as — A^, model 2 as — A'^, while 
model 3, —A. By GKN's general confinement criteria, 
model 1 and model 2 are super-confining, while model 3 
is critically confining. In the ultraviolet limit, the three 
models' dual gravity description have the same asymptot- 
ical background — AdS space-time plus running dilaton 
field. In the infrared region, model 1 and 2's dual space- 
time scale factor approaches to 1 as A — > cxd, but model 
3's scale factor shrinks to zero size as A — > oo. 

For future directions, we think looking for finite tem- 
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perature solutions to these models and the relevant phe- 
nomenological applications of them may be interes ting 
subjects, see e.g. references [s^l [s^l [Si [IHl [131 [23 
[3^. We hope to come back on this point in the near 
future. 



When a < 0, by definition, we can easily prove the fol- 
lowing iteration formula. 
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T[a, z] 



a<0 e z 



-ria + l,zl 



a(a -|- 1) a{a + 1) 



(A.4) 

r[a + 2,zl 



APPENDIX 



APPENDIX A: INCOMPLETE GAMMA 
FUNCTION 

This appendix collects some basic formulas of the in- 
complete Gamma function which is needed in the analyt- 
ical derivation of the short separation limit of the heavy 
quark potentials. The content is based on the help docu- 
ment of software Mathematica and some simple variable 
replacements and derivations. 

The incomplete Gamma function Fla, zl is defined as 



ria,zl 



(A.l) 



When a > 0, we have the following expansion formulae. 



a - 1 



ria, zi rial 



z I 

a 



1 + a 2(2 + a) 



z^ z 

(A.2) 

+ 0(z3)) 

(A.3) 



For integrations like \~Pe~'^dX, we can translate 
it into incomplete Gamma functions through variable re- 
placements 



-dX = 



X-Pe-^d\ = 



1 

cJ 
(A.5) 



'dX 



(A.6) 



C- 



ce 



[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2 (1998) 231- 
252, hep-th/9711200 

[2] S. S. Gubscr, I. R. Klebanov, and A. M. Polyakov, 
Phys. Lett. B428 (1998) 105-114, jhep-th/9802109 
E. W itten, Adv. Theor. Math. Phys. 2 (1998) 253-291, 
| hep- th/9802150; Adv. Theor. Math. Phys. 2 (1998) 505- 
532^ hep-th/9803131 

[3] I. R. Klebanov, Edward Witten, Nucl. Phys. B536 
(1998) 199-218, hep-th/9807080 ; I. R. Klebanov, N. A 
Nekrasov, Ncul. Phys. B574 (2000) 263,'hep-th/9911096j 
I. R. Klebanov , A. A. Tseytlin, Nucl. Phys. B578 (2000) 
123, hep-th/0002159; I. R. Klebanov, M. J. Strassler, 
JHEP 0008 (2000) 052, hep-th/0007191 

[4] L. Girardello, M. Petrini, M. Porrati and A. Zaffaroni, 
Nucl. Phys. B 569, 451 (2000), 9909047v2 [ hep-th] 
L. Girardello, M. Petrini, M. Porrati and A. Zaffa- 



roni, JHEP 9905, 026 (1999), |hep-th/9903026| L. Gi- 
rardello, M. Petrini, M. Porrati and A. Zaffaroni, JHEP 
9812, 022 (1998), 9810126 [hep-th] » D. Z. Freedman, 
S. S. Gubser, K. Pilch and N. P. Warner, Adv. Theor. 
Math. Phys. 3, 363 (1999), i9904017v3. D. Z. Freed- 
man, S. S. Gubser, K. Pilch and N. P. Warner, JHEP 
0007, 038 (2000), hep-th/9906194]j. Polchinski and 
M. J. Strassler, 0003136v2 [hep-thTl J. Babington, 
D. E. Crooks an d NT J. Ev ans, Phys. Rev. D 67, 066007 
(2003), 0210068 [hep-th] : N. R. Constable, R. C. My- 
ers, JHEP 9911,020 (1999), 9905081v4; J. Babington, 
D. E. Crooks and N. J. Evans, JHEP 0302, 024 (2003), 
0207076 [hep-th] 
15] T. Sakai, S. Sugimoto, Prog. Theor. Phys. 113 
(2005) 843-88 2, hep-th/0412141v5 Prog. Theor. Phys. 
114:1083-1118 |hep-th/0507073v4[ ~ 



9 



[6] 



A. Karch, E. Katz, JHEP 0206 (2002) 043, 
|hep-th/020 5236 D. Arean, D. E. Crooks, A. V. 
Ramallo, JHEP 0411 (2004) 035, hep-th/0408210; P. 
Ouyang, Nucl. Phys. B699 (2004) 207, hep-th/0311084 
T. S. Levi and P. Ouyang, hep-th/0506021 ; T. Sakai, J 
Sonnenschein, JHEP OSOQ (2003) 047, hep-th/0305049 
S. Kuperstein, JHEP 0503 (2005) 014, hep-th/0411097 

B. A. Burrington, J. T. Liu, L. A. Pando Zayas and D. 



[7] 



[8] 
[9] 

[10] 

[11] 

[12] 
[13] 

[14] 
[15] 



[16] 
[17] 

[18] 
[19] 
[20] 



Vaman, JHEP 0502 (2005) 022 |hep-th/0406207| B. A. 
Burrington, V. S. Kaplunovsky and J. Sonnenschein, 
'708.1234[hep-th]'; R. Apreda, J. Erdmenger, D. Lust 
and C. Sieg, JHEP 0701 (2007) 079 hep-th/0610276 ; C. 
Sieg, JHEP0708 (2007) 031 0704.3544[hep-th] 
R. Casero, C. Nunez and A. Parcdcs, Phys. Rev. D73: 
086005,2006, e-Print: hep-th/0602027 ; E. Caceres, R. 
Flauger, M. Ihl and T. Wrase, JHEP 0803:020,2008, e- 
Print: larXiv: 0711 .48781 F. Benini, F. Canoura, S. Cre- 
monesi, C. Nunez, A. V. Ramallo, JHEP 0709:109,2007, 
e-Print: arXiv : 0706 . 1238 1 

J. Erdmenger, N. Evans, I. Kirsch, E. Threlfall, 
10711.4467 [hep-th] , 

J. Polch inski, M. J. Strassle r, Phys. Rev. Lett. 88 (2002) 
031601, ■hep-th/0109174vl' 

S. J. Brodsky, G. F. dc Tcramond, Phys. Rev. Lett. 96 : 
201601 (2006), hep-ph/0602252 
A. Krikun, 0801 . 4215 [hep-th] 

L. D. Rold, A. Pomarol,"NucL Phys. B721 (2005) 79-97, 

hep-ph/0501218 

hep-th/0702205vl 

J. Erhch, E. Katz, D. T. Son, M. A. Stephanov, Phys. 
Rev. Lett. 95 (2005) 261602, hep-ph/0501128 
K. Ghoroku, N. Maru, M. Tachibana and 
M. Yahiro, Phys. Le tt. B 633, 602 (2006), 
I [arXiv:hep-ph/05 10334] | 

A. Karch, E. Katz, D. T. Son, M. A. Stephanov, Phys. 

Rev. D74: 015005, 2006, hep-ph/0602229v2 

O. Andrecv, V. I. Zakharov, Phys. Rev. D73 : 107901, 

2006 hep-t h/060 3170 ; 

Phys.Rev. D74l'200 6) 025023, |hep-ph/060420"4l 

E. Shuryak, 'hep-th/060 5219vl] 

F. Bigazzi, R. Casero, A. L. Cotrone, E. Kiritsis, A. Pare- 
des, JHEP 0510 (2005) 012, hep-th/0505140j 

U. Giirsoy and E. Kiritsis, eprint: 0707.1324. 
U. Giirsoy, E. Kiritsis and F. Nitti eprint: 0707.1349. 
P. Boucaud, F. De Soto, A. Le Yaouanc, J. P. Leroy, J. 
Micheli, H. Moutarde, O. Pene, J. Rodriguez-Quintero, 
,JHEP 0304(2003)005, |hep-th/02 12192] 



[21] 
[22] 
[23] 

[24] 

[25] 



[26] 



[27] 



[29] 

[30] 
[31] 



M. Frasca, e-Print "arXiv : 0802.1183 

T. A. Ryttov, F. Sannino, arXiv: 0711.3745 1 

G. S. Bah, P hys. Rept. 343 : 1-136 (2001), 

'hep-ph/0001312 

E. Eichten, K. Gottfried, T. Konoshita, K. D. Lane and 
T. -M. Yan, Phys. Rev. D17 (1978) 3090, D21 (1980) 203. 
J. M. Maldacen a, Phys. Rev. Lett. 80, 4859 (1998), 
|hep-th/9803002j 

S. J. R ey and J. T. Yee, Eur. Phys. J. C22, 379 (2001), 
|hep-th /9803001v3 

Y. Kinar, E. Schreiber, J. Sonnenschein, Nucl. Phys. 
B566 (2000) 103-125, 9811192; 

Y. Kinar, E. Schreiber, J. Sonnenschein, N. Weiss, Nucl. 
Phys. B583 : 76-104(2000) , hep-th/99 1 1 123v2 , 
D. Bak, A. Karch' and L. O* YafTe, 
le-Print: arXiv: 0705. 09941 

C. D. White, Phys. Lett. B652 : 79-85, 2007, 
e-Print: hep-ph/0701157 

J. P. Shock, F. W u, Y. -L. Wu, Zh. -F. Xie, JHEP 
0703:064,2007. |e-Print: hep- ph/0611227 



S. Weinberg, Cambridge University Press, 1996. 
G. S. Bah and K. Schilhng, Phys. Rev. D47: 661- 
672, 1993, e-Print: hep-lat/9208028 ; G. S. Bah, 
K. Schilhng and A. Wachtcr, Phys. Rev. D56: 2566- 
2589, 1997, e-Print: hep-lat/9703019 ; S. Necco 
and R. Sommer, Nucl. Phy s. B622: 328-346, 2002, 
|e-Print : hep-lat/0108008| 
[32] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Rev. 
Mod. Phys. 77: 1423,2005, e-Print: hep-ph/0410047; 
N. Brambilla, J. Ghiglicri, A. Vairo and P. Petreczky, 
le-Print : ~ arXiv : 0804 . 0993 By Quarkonium Working 
Group (N. BrambiUa ct al.). Published as CERN Yellow 
Report, ^e-Print: hep-ph/0412158,, 
[33] C. P. Herzog, Phys. Rev. Lett. 98:091601,2007. 

"hep-th/0608151' 

[34] Phys. Lett. B645 (2007) 437-441, | hep-ph/0607026| 

[35] E. Nakano, S. Tcraguchi, W.-y Wen, Phys. Rev. D75: 

085016,2007. hep- ph/060 8274 
[36] U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, eprint 

arXiv: 0804. 08991 
[37] S. S. Gubser, A. Nellore, S. S. Pufu and F. D. Rocha, 

e-Print: arXivT08643950 
[38] A. Clicrman, T. D. Colicn and E. S. Werbos, e-Print: 

arXiv: 0804. 1096 
[39] N. Evans and E. Threlfall, e-Print: I arXiv: 0805. 09561 



